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Abstract 

We derive Fourier's law for a completely coherent quasi one-dimensional 
chaotic quantum system coupled locally to two heat baths at different 
temperatures. We solve the master equation to first order in the temper- 
ature difference. We show that the heat conductance can be expressed 
as a thermodynamic equilibrium coefficient taken at some intermediate 
temperature. We use that expression to show that for temperatures large 
compared to the mean level spacing of the system, the heat conductance 
is inversely proportional to the level density and, thus, inversely propor- 
tional to the length of the system. 

PACS numbers:44.10.-|-l, 66.25. +g 

1 Introduction 

In classical physics, Fourier's law states that the heat conductance is inversely 
proportional to the length L of a physical system. For a system S coupled at 
either end to two reservoirs at different temperatures Ti and T2 with T2 > Ti, 
the heat conductance C is defined by writing the heat current through the 
system as C(T2 — Ti). The law C cx is intuitively obvious when one thinks 
of as the resistance of a macroscopic system against heat (i.e., phonon) 
transport. The system can be thought of as consisting of building blocks each 



with its own resistivity against heat transport. The resistance is the sum of 
these resistivities and grows hnearly with L, resulting in C cx 1/L. 

Does that law also apply when S* is a quantum system? The argument just 
given suggests that the answer depends on the degree of coherence of the system. 
Consider, for instance, a system consisting of K building blocks of length I each 
so that L = Kl, and assume that at the interface between neighboring blocks 
the system is coupled to the outside world so that quantum coherence between 
neighboring blocks is destroyed. Then the classical argument given above applies 
and C is inversely proportional to L. 

Here we are interested in the heat conductance of a fully coherent quasi one- 
dimensional quantum system S. No sources of decoherence are present other 
than the decoherence due to the coupling of S to the reservoirs at either of its 
ends. Phonons travel coherently from one end of S to the other. Because of 
that coherence, the inverse heat conductance cannot be calculated by adding 
the resistivities of individual blocks in the manner described above, and it is 
not clear why C should be proportional to L~^. The example of electron trans- 
port through mesoscopic samples coupled to two external leads actually suggests 
independence of C on L. That is seen as follows. At sufficiently low tempera- 
ture, electron transport is fully coherent [T. The electrical conductance cannot 
be calculated by adding the resistivities of parts of the system. Rather, for 
non-interacting electrons the conductance as given by the Landauer-Biittiker 
formula is the sum of squares of elements of the scattering matrix, i.e., of the 
quantum-mechanical transition amplitudes connecting the two leads. That for- 
mula embodies full quantum coherence of the mesoscopic system. The elements 
of the scattering matrix do not display a systematic dependence on the length 
of the system. Thus, the analogy with electron transport suggests that in a fully 
coherent system, C should be independent of L. On the other hand, numerical 
studies of several small quantum systems (typically spin chains) coupled to two 
reservoirs strongly suggest that C is indeed inversely proportional to the length 
of the system [Idll]. 

In the present paper we offer an analytical resolution of the resulting para- 
dox. We focus attention on chaotic quantum systems. We thereby exclude 
both integrable systems and disordered systems with localization. We do so 
because it is known that some integrable systems do not comply with Fourier's 
law. Moreover, it is clear that in disordered systems with localization, the heat 
conductance decreases exponentially with length. We later identify those of our 
assumptions that fail for localized systems. Starting from the master equation 
which describes the coupling of a quantum system to two reservoirs (heat baths) 
at different temperatures, we use a perturbation expansion in the temperature 
difference to derive in first order an expression for the heat current. We show 
that the heat conductance C in that expression can be written as a thermody- 
namic equilibrium coefficient calculated at some intermediate temperature Tq 
even when the conducting system itself is not in thermodynamic equilibrium. 
We show for a wide class of quantum systems that except for a normalization 
factor, the resulting expression for C is indeed independent of L as suggested 
by the analogy with mesoscopic electron transport. The normalization factor 
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has the form 1/ X^m 6xp[— £'m/(fcBTo)] typical for equhbrium systems. Here 
are the eigenenergies of the system S and fcs is the Boltzmann constant. It is 
easily seen that the sum over m is linear in the total level density of S and, 
thus, grows linearly with L. For C that yields an inverse dependence on L in 
agreement with the results reported in Refs. 

Work on the heat conductance in quantum mechanics has a long history, 
starting with Refs. O [6l [7], and many papers thereafter. The length depen- 
dence of the heat conductance C and the influence of quantum coherence on 
the value of C were discussed in none of these early works, however. In the 
related electrical-conductance problem, the importance of quantum coherence 
was recognized only in the 1980s jT|. During the last decade Fourier's law in 
quantum mechanics has been intensely studied, especially for quantum spin 
chains. In addition to the papers cited above we mention Refs. [8j |9l [TOl [11] 
and the review [3] where further references may be found. The present paper is 
based on a random-matrix model. The results derived within that model and 
obtained as ensemble averages are generically valid for quasi one-dimensional 
chaotic quantum systems. 



We consider two heat baths (reservoirs) labeled 1 and 2 with temperatures Ti 
and T2 > Ti coupled to a quantum system S. The Hamiltonian Hs of the system 
has eigenfunctions |m) and eigenvalues Em- The occupation probabilities of the 
states \m) are denoted by P^- The stationary state of the system S is described 
by the master equation 



The derivation and conditions of validity of Eq. ([T]) have been thoroughly dis- 
cussed in the literature, see for instance Ref. [12]. To fix our notation, we give 
a brief derivation of Eq. ([T]) in the Appendix. We will show that the Pm are 
uniquely determined and, for Ti ^ T2, in general differ from an equilibrium 
distribution. We will use Eq. ([T]) to calculate the stationary heat current and, 
from that, Fourier's law and an expression for the heat conductance. In doing 
so we exclude systems and/or couplings with pathological properties. Examples 
would be a system consisting of two uncoupled parts, or a localized system, or 
a system with coupling matrices Wmn where at least one is close to diagonal in 
the eigenvalue representation of Hs- 

It is clear from the outset that within the framework of Eq. ([T]) a temper- 
ature gradient within the system 5" cannot exist since the temperature of the 
system S is defined in terms of the occupation probabilities Pm and these are the 
same throughout the system. That fact reflects quantum coherence throughout 
the system S- This statement does not preclude the possibility that a suitably 
defined local energy expectation value possesses a non-zero gradient, however. 



2 Master Equation 



m 
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For such a situation to arise, the P„i rnust differ from the thermodynamic equi- 
hbrium values in Eq. (|4]) below. 

Explicit values of the coupling matrix elements Wmn{Ti) with i = 1,2 are 
worked out in the Appendix. According to Eq. (1551) these obey 

W^i (TO = X^l exp{ (A/2) {E„, ~ E^)} (2) 

where Pi = l/T,;. We put the Boltzmann constant equal to unity. Here Xnm is 
real, independent of temperature and, in contrast to Wmn, symmetric, Xnm = 
Xl^. Eq. © implies 

WiiiiTi) = exp[A(i?™ - E^i^{T{) . (3) 

Eq. ([3]) immediately yields the form of the normalized solutions P„ of Eq. ([T]) 
for the equilibrium case, Ti = T2 = T. With (3 — l/T these are given by 

^ exp[-pEn] 
" E™exp[-/3£;,„] • ^ ' 

The equilibrium distribution Q is independent of the values of the elements 
of the coupling matrices X'-'^\ Can such independence also be expected for the 
non-equilibrium case Ti ^ T2? For simplicity we take = A:(2) ^ X and 
note that the matrix X is symmetric. Eq. ([T]) takes the form 

^X„„| (^exp[(/3i/2)(£;,„ - En)]+exp[{f32/2){E,n - En)]^P,n 

exp[(/?i/2)(^„ - En,)] + exp[(/32/2)(S„ - E„,)]^Pn^ = . (5) 

For this equation to hold independently of the values of the Xnm, each of the 
coefficients multiplying Xmn must, in general, vanish individually. It is seen 
immediately that that is possible only for Ti = TJj- The developments in Sec- 
tions 13.11 and 13.21 below show that the argument is somewhat simplistic and 
must be refined. Nevertheless we conclude that in the non-equilibrium case the 
values of the occupation probabilities Pm depend, in general, on the values of 
the X^'^ 

To determine the stationary heat current / through the system S we muhiply 
Eq. ([T]) with En and sum over n. We obtain 



mn 



EnWill{T,)Pn, - J2 EnW^liT^)Pn 

mn 

= -J2 EnWi'2iT2)Pm + J2 EnWiSAT2)Pn ■ (6) 
mn mn 

Eq. ^ expresses energy conservation. We interpret the negative left- and the 
positive right-hand side of Eq. ^ as the heat current / (energy per unit time 
transferred to the system from bath 2 or from the system to bath 1, respectively). 
That interpretation is in keeping with the fact that the master equation ([T]) is 
the stationary form of a more general equation for the time derivative of the 
occupation probabilities Pm of the states m. 
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3 Perturbation Expansion 



According to Fourier's law, the heat current / is proportional to the tempera- 
ture difference between the two heat baths. This suggests using a perturbative 
approach in powers of 5T ~ (1/2)(T2 — Ti) > 0. We accordingly expand the 
quantities appearing in Eq. around some intermediate temperature Tq that 
obeys Ti < To < T2 in powers of 5T up to and including terms linear in 5T . 
The choice of Tq plays an important role in the calculation. 

Expanding the coefficients W^'^iT,) at To in powers of 5T = (1/2)(T2 - Ti), 
we use 

^=.^(£;„-K,„)W^W , (7) 

see Eq. ([2|) . We likewise expand the solutions Pm of Eq. ([1]) around the equilib- 
rium solution Q at T = To, omitting the normalization factor, 

P„ = exp[-/3oS™](l + 5P,n) (8) 

where /3o = I/Tq. By definition of the equilibrium solution, the terms of zeroth 
order in 5T in Eq. ([T]) mutually cancel. For i = 1, 2 we define the real symmetric 
matrices 

B^± = exp[-(/?o/2)(i?„ + i?„)]X«„ . (9) 



and the vectors 



We note that 



^0 



^A«=0. (11) 

m 

The master equation takes the form 

A A tin ^ <i m ' 



(12) 



This is a set of inhomogeneous linear equations for the unknown quantities 
6Pn- The homogeneous equations possess the non-trivial (equilibrium) solution 
5Pn = 1 (all n). According to Eq. pT|) the inhomogeneity is orthogonal upon 
that solution. Therefore, the inhomogeneous equations possess a unique solution 
5Pm with 

^ SP,n = . (13) 



E 



That is the solution we study in the sequel. 

We consider consecutively three cases: (i) The coupling matrix elements are 
equal, Xmn = Xmli] (ii) the coupling matrix elements are similar, Xmn = aXmn 
with a > 0; (iii) the coupling matrix elements are dissimilar so that neither case 
(i) nor case (ii) applies. 
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3.1 Symmetric Coupling 

For the symmetric case with equal couplings, Xmn — Xmn — Xmm we choose 
the intermediate temperature 

To = i(Ti+r2). (14) 

We have Bmn — Bmn — Brnn SO that J2ii~y — ^ ^^''^ Therefore, 
the solution of Eq. vanishes identically, and the solution of the master 
equation ([T]) is given by the equilibrium solution Q with 1 //3 taken at T = Tq 
even in the non-equilibrium case Ti ^ T2, up to and including linear terms 
in (5T. The system S is in thermal equilibrium at the mean tempertaure Tq. 
There cannot exist a gradient in the local expectation value of the energy. A 
non-vanishing non-equilibrium solution of the master equation ([l} must be of 
second order in bT . 

An explicit expression for the heat current is obtained by using for P,„ the 
equilibrium solution (including the normalization factor) and by expanding the 
left-hand side of Eq. ^ up to terms linear in JT. That gives 

^ = E(^™ - ^")'^™„) . (15) 

This, in a very general form, is Fourier's law. The heat conductance is given by 
an equilibrium property of the system: The coefficient J2mni-^m ~ EnYBmn is 
calculated at the equilibrium temperature Tq. The choice (fT4| of To is justified 
by hindsight: It is easy to check that for a different choice, we do not find 
5Pm = for all m. That shows the importance of the choice of Tq: Apparent 
deviations from thermal equilibrium indicated by nonvanishing 5PmS may only 
be caused by an improper choice of Tq. 

3.2 Similar Couplings 

We turn to case (ii) where Xmli = aXmn with a > 0. Again, the choice of Tq 
is important to minimize apparent (but unreal) deviations of the system from 
thermal equilibrium. Thus we write 

To = aTi + (1 - a)T2 (16) 

with < a < 1 and determine a from the condition that the solutions SPm of 
the linearized master equation (jl2p vanish. That equation now takes the form 



[-2(l-a)A«+2aAf ]5T + ^^S« JP,„- ( ^^bW Up„ = . (17) 

Because of the definition pO|) and the condition Xmn = aXmn the inhomogene- 
ity vanishes if we choose a = a/{l + a) and, thus. 

To = {aTi + T2) . (18) 
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Thus, the system S is in thermal equihbrium at the temperature Tq given by 
Eq. (jlSp . The equihbrium temperature is shifted away from the arithmetic mean 
of Ti and T2 toward the temperature of the heat bath with the stronger couphng 
to the system. That is physically very plausible. The heat current is given by 

= J^^(i^,„_i5„)2s(2A ^r. (19) 

The same conclusions as in Section [3 . 1 1 applv : The system is in thermal equilib- 
rium at temperatute Tq. The heat conductance is evaluated at the equilibrium 
temperature. A gradient of the local energy expectation value and deviations 
from equilibrium must be of second order in ST. 



3.3 Dissimilar Couplings 

In the case of dissimilar couplings we expect that the difference will be rather 
small and statistical in nature. Indeed, the derivation of the master equation in 
the Appendix is based upon the assumption that the operators Q^*^ that couple 
the system S to either heat bath are local operators that act on either end of 
S. In describing such coupling operators explicitly, it is commonly assumed 
that they are made up of local position and momentum operators. For chaotic 
systems we expect the matrices X'^^ and X^'^^ to be similar. It is not difficult, 
however, to imagine systems where the matrices AT*-^' and A'^^ differ strongly. 
That is the case, for instance, for localized systems. Then A*^^' couples most 
strongly to states localized near one end of the sample. Such states do not 
couple significantly to A^^^ at the other end, and conversely for X'^^K A drastic 
and systematic difference between the two coupling matrices invalidates our 
treatment. 

The case of dissimilar couplings is more complicated than the cases of equal 
and similar couplings and requires some algebra. It turns out that the occu- 
pation probabilities P„i differ from an equilibrium distribution. However, it is 
possibe to define a temperature Tq such that the heat conductance is given by 
an equilibrium expression calculated at Tq. 

We use the ansatz Eq. and obtain Eq. (fT7|) . But now the inhomogeneity 
in Eq. I\17^ does not vanish for any value of a, and it is necessary to determine 
the solutions 5Pm oc 5T of that equation. We define Bnm = J2i Bum and the 
real symmetric matrix i3„,„ = Bnm — Snm Bnk- Eq. (|17p takes the form 

[-2(1 - a)A';p + 2aA^^^6T + ^ B^mSPm - . (20) 

We have J2m Bnm = 0. Therefore, the matrix Bnm possesses one vanishing 
eigenvalue, Ai say, with associated eigenvector {1, 1, . . . , 1}-^. The matrix B 
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can be diagonalized by a real orthogonal matrix O, 

[OBO^]„m ~ Smn\n ■ (21) 

The eigenvector {1,1,..., I}-'" occupies the first column of the matrix , and 
the vector {1,1,...,1} occupies the first row of the matrix O. We assume 
that all eigenvalues Am with m > 2 differ from zero. We multiply Eq. (|20p 
from the left with the matrix O, define the vectors Am = J2k ^mfc^fe ^''^d. 
SP^ = J2k OmkSPk, and observe that A^*' = (see Eq. ^) and that 5~Pi = 
(see Eq. dfe])). For m > 2 that yields 

[-2(1 - a)i(^) + 2aA\l^]5T = X^SP^ . (22) 

Solving for 5Pm and transforming back to 5Pm we find 

5Pra = StY,(T. (^km^Okn) [-2(1 - a)4'^ + ^aA^n^] . (23) 

In order to use that result in Eq. we normalize the solutions defined by 
Eq. (|S]) so that J2m Pn — 1 and expand the resulting expression in powers of 
SPm, keeping only terms up to first order. For P,„ that yields 

eM-PoE,n] A , ,p Efcexp[-/?ogfc]JPA 
" E„ exp[-/3oP„] " exp[-/3o£;d j " ^'^^ 

Insertion of Eqs. into ([M)) gives the occupation probabilities of the states m. 
The solutions 5Pm do not vanish identically for any choice of a, and there is no 
choice of temperature Tq for which the system would be in thermal equilibrium. 
Thus, the local expectation value of the energy may not be the same throughout 
the system, and it may possess a non-zero gradient. We have not investiated 
that possibility. 

We turn to the heat current / defined in Eq. In terms of the solu- 

tions ([23| / is given by 

^ = ^^-FTfS2(l-a)^(£;,„-P„)2i?W 

T,n exp[-/3o-B„] V ^ 



E„exp[-/3oii;„] V2r2 

-J2{Em - E.^)B^^lSP,A . (25) 



To simplify these expressions we consider the last term in the first of Eqs. 
With the help of Eqs. ([TT|) and ((23|) that term can be written as 



mn ^ k>2 k>2 



Y^iEm - E^)B^^JP^ = 2TS6T ( a ^ ^if - (1 " «) E 4^' Yk^i'' 



(26) 
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For equal couplings we have A^j^^ = A^j^^ and the right-hand side of Eq. (|26l) 
vanishes for a — 1/2. For similar couplings we have A'j^^ = aA^^'^ and the 
right-hand side of Eq. ([26]) vanishes for a = 1/(1 + a). In the present case of 
dissimilar couplings, the expression in big round brackets is linear in a and, 
therefore, vanishes at some uniquely defined value ai. That value defines via 
Eq. (HH) a temperature T^^\ We expect that Tfj^^ obeys Ti < T^^^ < T2 and 
that, therefore, ai obeys < ai < 1. Values of ai outside the interval [0, 1] 
would be physically implausible. Similar considerations apply to the second of 
Eqs. (pSj) . Again there exists a value of a2 and an associated temperature Tg ' 
for which the term linear in SPn vanishes. We cannot prove, however, that 
a2 = ai, or that T^^' = T^'^'. 

Choosing Tq = T^'^^ in the first and Tq = T^^"^^ in the second of Eqs. ^ we 
obtain for the heat current 



/ = 



E„exp[-/3«i?„]V(r(i))2 

1 ^ I 



In these equations, the coefficients B*-*-* are taken at the temperatures Tq*' with 
i = 1,2, and = 1/Tii'^- We cannot show that T^'^^ = T^'^'^ although that 
would be physically most plausible. The heat conductance in both expres- 
sions (l?fl) is evaluated at some equilibrium temperature Tq*'' even though the 
system is not in equilibrium. 

4 Length Dependence of the Heat Conductance 

In Section [3] we have shown that in Fourier's law, 

I 6T , (28) 

the heat conductance C can always be written as an equilibrium coefficient of 
the form 

^ ^ ^ TTTITT - ^")' cxp[-(/3o/2)(i?™ + £;„)]X„„ . (29) 

That expression may be viewed as a special case of the Green-Kubo formula O 
[6]. Here 7 is a numerical coefficient of order unity, and Tq is a suitably defined 
temperature. The symmetric matrix Xmn describes the coupling of the quantum 
system S to one of the two heat baths. We now investigate how the heat 
conductance C depends on the length L of the system S. We consider a linear 
chain or a piece of wire or some other quasilinear system, all of length L. We 
model the system in terms of random matrices, thereby assuming that it is 
chaotic. We show that for any such device, coupled at either end to heat baths 



\2 



a^Y^iE^ - E^fB(^l{Tl,'^))ST . (27) 
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with different temperatures, the heat conductance is inversely proportional to 
L as is the case in classical physics, and as is found to be the case for small spin 
chains [Hll]. 

To this end we rewrite the terms appearing in Eq. (j29l) . using the defini- 
tion ®, 

^exp[-/3oS™] = I AE eM-PoE]Y,5{E-Em) , 

m m 
^(-Em — EnfBmn — ^(i?m — E^f exp[— /3o(£^m + En)]Xmn 



J dei y"de2(ei -e2)^exp[-/?o(£i +62)] 
x27rAo cxp[(ei - £2) 



mn 



<5(£i - Er^)S{e2 - Er,) . (30) 



We have used Eqs. ^ and ([55)1 and written the expressions appearing in the 
heat conductance C in terms of length-independent energy integrals, in terms 
of IQmnP, and in terms oi 5{E — Em)- The constant ^0 measures the strength 
of the couphng of the system to the heat bath, see Eq. (|46)) . By definition the 
operator Q couples the surface of the system S locally to the heat bath and does 
not depend on L. Any length dependence of the terms in Eq. pop is due to sums 
involving 5[E — Em)- Such sums appear differently in the first and in the second 
of Eqs. ([50]) . In the first, the level density 5{E — E^) itself appears as an 
independent factor while in the second, the sums involve the squared matrix 
element |QmnP of Q. In condensed-matter physics, expressions of similar form 
are referred to as the local density of states and in nuclear physics, as the 
strength function of the operator Q. It is obvious and confirmed below that the 
density of states appearing in the first of Eqs. ([30]) increases linearly with the 
length L of the system S. It remains to show that the length dependence of the 
last term in the second of Eqs. ((30)l is negligible. We do so for a large class of 
quasi one-dimensional systems S. 

We exploit the fact that Q is local, i.e., it acts only on one surface of S. 
Let P denote an orthogonal projector obeying P ^ P^ = P^ that projects onto 
a suitably chosen set of states on that surface such that Qmn = {m\Q\n) = 
{mP\Q\Pn). We investigate the length dependence of the expression 

y^\Pm)S{E-Em){mP\^--P^—^-—P (31) 

m 

where E~^ carries an infinitesimal positive imaginary increment and where Hs is 
the Hamiltonian of S. Expression ([51]) differs from the expression for the total 
level density 

p{E) = E ^(£; - Em) - -^5 Trace-^^^-^ . (32) 
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We show that the projector P in Eq. pip has a profound influence on the length 
dependence. 

To model the length dependence of S we think of S as consisting of K blocks 
of fixed length / each, labelled by a running index k — 1, . . . ,K. The length 
of S is then given by if = Lk, and the dependence on length is converted into 
a dependence on K. For Hs we use a matrix representation. Let i?*^*^^ be 
the Hamiltonian matrix for block k. Only neighboring blocks are coupled by 
Hamiltonian matrices lyC^'^+i) , and we have 



Hs = 



( iJd) 




V ... 



i/(2) 
M/(23) 







1^(23) 
i7(3) 
Vl/(34) 






1^(34) 
i/(4) 



(33) 



where the upper indices range from 1 to fc. 

Let the block with k = \ carry the surface on which Q acts. Then we have 
PHsP = PH^^^P. In the first block, we introduce a complete set of orthonormal 
basis states where /i = 1, 2, . . .. The basis is chosen such that only the first 
n states are surface states so that P\^i) = for fj, > n. Then P ~ X]J^=i 
Instead of the strength function in Eq. (|3T|) we consider 



1(11 / 

TT ^ ' E+ - Hs 



ID 



(34) 



Without loss of generality we confine ourselves here to the diagonal element of 
the propagator with respect to the state with fj, = 1. That choice is arbitrary, 
any other value of fi with fJ. < n would give the same result. We consider 
only diagonal elements because the non-diagonal elements vanish on average 
as a consequence of the statistical assumptions introduced below. With Ei 
(l|if(i)|l) and K 
written as 



h[]^ for 1^ > 2, the matrix in Eq. ([33]) is explicitly 



El 



(35) 



Here v, p> 2. We put Ei ~ Q and comment on that choice below. 

Analytical progress is possible upon introducing statistical assumptions on 
the matrix elements of Hs- We assume that the matrices H'^^^ with k = 
1, 2, . . . , fc are uncorrelated and are members each of the Gaussian Orthogonal 
Ensemble (GOE) of random matrices of dimension N ^ with the limit N ^ oo 
eventually taken ^3j. Ensemble averages are denoted by angular brackets. The 
matrix elements are Gaussian-distributed random variables with mean values 
zero and for k,l,m = 1, 2, . . . , if obey 



\2 

(1 - SikjSklJ^i^tiv'^un' + 5^ti'Svy') 



w 

'n 



{Skk'Sw + 5kv5ik'){5^^'6ppi +5pp>5fj_>p) , 



X-^-^pa p'a' I 



(36) 
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Because of these assumptions the average spectrum of each of the matrices 
ijC^) with k > 2 has for A'' — >■ oo the shape of a semicircle with radius 2A. 
The coupling between neighboring blocks mediated by the matrices W is char- 
acterized by the strength parameter w. It is physically obvious that w < X. 
For the first block we introduce a special notation. We assume that the first 
of Eqs. applies in form also to the case where k ^ I ~ 1 but does so only 
for ii.v^n'.y' > 2. We assume that the in Eq. (|35|) are Gaussian random 
variables with mean value zero, not correlated with the other matrix elements, 
and obey {V^Vi,) — v'^S^i,. We do so in order to display explicitly the role of the 
coupling of the surface state to the rest of the system, without any restrictions 
on the value of w^. For N ^ oo the spectrum of H^^^^ is also of semicircular 
form but the surface state plays a distinct role. 

The random-matrix model introduced in Eqs. (1361) describes the generic 
properties of chaotic quasi one-dimensional quantum systems. It is the most 
general model we can think of to describe such systems. The same model has 
been widely used to describe electron transport through disordered mesoscopic 
samples [T¥l ITS] . 

We work out the ensemble averages of expressions (|32]) and (p4)) in the frame- 
work of the random-matrix model defined in Eqs. p6p. To this end we calculate 
the average Green function G{E) = {{E~^ — Hs)~^) of the system. Because of 
the Gaussian distribution of the elements of Hs, that function obeys for — 7> cx) 
the Pastur equation [TB] 

EG{E) - 1 + {HsG{E)Hs)G{E) . (37) 

In view of our statistical assumptions, all non-diagonal elements of G{E) (with 
respect to both block index k and running index /i) vanish. We write G^^l{E) 
for the diagonal elements in block k and define G^'^^ — (l/N) G^^\ In block 
1 the sum runs from /i = 2. We consider the matrix element G^^^ (E) separately 
because the expression in Eq. (l34l) is given by {—l/Tr)'^G^ii{E). According to 
Eq. dSIl) we have for A^ > 1 and fc = 1, 2, . . . , AT and with G(°) = = G'^+i) 

EG[^j'{E) = 1 + v^NG^g[^j' , 

EG^''\E) = 1 + [A^gC^) -I- iw^Gf^-i' + w^G'-'^+^^jG'^''^ . (38) 
Therefore, g']^-^ (E) is given by 

We observe that the factor Nv'^ represents the total strength of the coupling 
of the surface state |1) to the system S. Since the state |1) couples only to 
states in block 1, that strength is independent of the length of S. The complex 
propagator G^^^ is obtained by solving the second set of Eqs. (p8l) . 

We solve these equations approximately. For A' ^ 1 and k somewhere in 
the middle of the range [1, AT], the form of Hs in Eq. (pS)) suggests that G^*^-* 
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changes slowly with k. We accordingly put G^'^^ = G'^'' = G^''^^\ With 



(A')' = A2 + 2w^ (40) 
the resulting quadratic equation for G^'^^ yields 




yG^'HE)^^-^,|l-{§^'^\ (41) 

The average spectrum in block k, proportional to the imaginary part of G^'^^ , 
retains the form of the semicircle. However, the range 4A' of the spectrum 
is increased compared to the value 4A that applies without coupling to the 
neighboring blocks {W = 0). The increase is independent of K, i.e., of the 
length of the system. 

For blocks near the either end of the system, i.e., A:- values close to 1 or 
K, the form of Hs in Eq. p3|) and the form of Eqs. ([38)) suggest that the 
solutions retain the form (|4T|) but with values of A' that are smaller than given 
in Eq. (|40|) . That statement is supported by taking K — 2 in which case we 
find (A')^ = A^ + w"^. We conclude that the solutions of Eqs. (pSj) yield average 
spectra of approximately semicircular shape with ranges A' that lie between 
VA^ + and VA^ + 2w'^ and that are independent of K. An estimate of 
the range AE of the spectrum of Hs confirms that conclusion. We use that 
AE w ^/{l/{KN))TTace{Hs)^. For > if > 1 we obtain AE « VA^ + 2w^, 
a result that is independent of K and consistent with the values for A' just 
mentioned. 

Using the form dH]) for G'-^^E) in expression p4l) and putting in G'^-'(i?) 
the energy argument equal to zero for simplicity we obain 

1 1 

(42) 



TT E + iv^N/X' ' 

The width 2v'^N/X of the strength function of the state |1) is essentially given 
by the total strength Nv'^ of the coupling of that state to the system S divided 
by the range of the spectrum of S. That is a physically very plausible result. 
Neither Nv'^ nor A' depend on the length of the system. The result (|42l) can 
easily be extended by taking into account the full form of the propagator in 
Eq. (|4T|) and by dropping the assumption Ei = 0. Instead of the form (|42]) one 
obtains an expression of Breit-Wigner form centered at Ei with a width and a 
level shift due to the imaginary and the real parts of G'^^\ respectively. That 
expression is also independent of K. The same statement obviously applies to 
the full expression (PT|) . 

We compare this result with the total level density of S given by Eq. (|32|) . 
From the definition of G'^'^' it follows that 

(p(£;))=--3^G«(£;) (43) 

k 
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and, with G'^^\E) almost independent of k, {p{E)) cc K cc L. We conclude 
that the double sum in the numerator of Eq. (|29l) is independent of L while the 
single sum in the denominator is linear in L. As a result the heat conductance 
C is inversely proportional to L. 

5 Summary and Conclusions 

For a completely coherent chaotic quasi one-dimensional quantum system cou- 
pled to two heat baths at different temperatures, we have solved the master 
equation up to first order in the temperature difference and obtained Fourier's 
law. For equal or similar couplings at both ends of the system we have shown 
that the heat conductance C in Fourier's law can always be written as an equi- 
librium coefficient. In the case of dissimilar couplings the same statement holds 
generically for chaotic systems but not, for instance, for localized systems. 

We have used that result to investigate the dependence of C on the length L 
of the quantum system. Intuitive arguments based on quantum coherence and 
the analogy with electron transport through mesoscopic systems both suggest 
that C be independent of L, in contrast to numerical evidence [U |3] showing 
that C oc L~^. We have resolved that discrepancy by showing that aside from 
a normalization factor, C is indeed independent of L. The entire length depen- 
dence of C is found to be due to the normalization factor and determined by 
the density of states. The latter increases linearly with L and yields C oc L~^. 

The length dependence of the remaining term in C is determined by that of 
the last factor in Eq. ([50)1 . That factor bears a close analogy to the spreading 
width in nuclear physics and to the local density of states in condensed-matter 
physics. In both cases, it is known that the values are not affected when the 
dimension of the system is increased. In the present case, we have used a 
random-matrix approach to model the length dependence. That approach yields 
generic results |13) (exceptions have measure zero with respect to the probability 
measure that defines the random-matrix ensemble). We have also used the 
essential fact that the coupling to the heat baths is local and linked to the 
surface. We have shown that that fact guarantees the length-independence of 
the relevant term in Eq. (l30t . 

Having resolved the discrepancy, we may turn the question around and ask: 
Why does electron transport in mesoscopic systems not show a similar length 
dependence due to a normalization factor? The answer lies in the tempera- 
tures at which heat transport and electron transport are considered. To insure 
quantum coherence, electron transport is experimentally studied close to zero 
Kelvin. At such low temperature, only the lowest eigenvalues Em of the system 
would contribute to C. The random-matrix model we use is viable only when 
the density of states is sufficiently high, i.e., when the temperature Tq at which 
C is evaluated, is very much larger than the mean level spacing near the ground 
state. The random-matrix model is not expected to account correctly for effects 
near the end points of the spectrum. 

Although C is inversely proportional to L in both classical and quantum 
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physics, the causes for that dependence are seen to be strikingly different. In 
classical physics and for systems of macroscopic size, lack of quantum coher- 
ence causes the total resistance to be the sum of the resistivities of subsystems 
and, thus, C oc L~^. For a fully coherent quantum system, the inverse length 
dependence of C is due to the linear increase with L of the level density of the 
system. 
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Appendix: Derivation of the Master Equation 

In Ref. [T7], the master equation for a quantum system coupled to a single 
heat bath was derived. Here we use the same derivation for a quantum system 
coupled to two heat baths with temperatures Ti and T2 to obtain Eq. ((TJ. We 
only sketch those parts of the derivation that differ from Ref. [17] and do not 
repeat here the discussion of the conditions under which the derivation holds. 
These are the same as in Ref. jl7| . 

The quantum system S has Hamiltonian Hs, eigenvalues Em, and eigen- 
functions |m), with m — 1,2,.... With i = 1,2 the two baths have Hamiltoni- 
ans Hi, eigenvalues ei* '' and eigenfunctions \ia), respectively, with a = 1, 2, . . .. 
The coupling between each of the baths i = 1,2 and the system has the form 
Wiz) ^ Q(^)]/(i). Here Q(*) are local operators that act on either end of the sys- 
tem S while the V^l' are uncorrelated random matrices to be specified below. 
The total Hamiltonian is 

i 

= Hq + W . (44) 

To define the ensemble we consider case (II) of Ref. [T^ and assume that the 
are Gaussian-distributed random variables with zero mean values and 
second moments given by 

{mia I | nib) (pic \ | qid) = [SmpSnqSacSbd + SmqSnpSadSbc] 

(laTO|T/F(i) \lbn) {2cp\ W^^^2dq) ^ . (45) 

The overbar denotes the average over the ensemble. The second of Eqs. (|45p 
shows that VF'-^' and ly(^) are uncorrelated. In the first of Eqs. (|45p we assume 
standard random-matrix properties. We assume that the coupling of the system 
S to the heat bath described by the matrices Q*-*' acts only on the (left or right) 
surface of the system and not on its volume. For the matrices y*^*' we assume 
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as in Ref. [H] that with z = 1, 2 



[V^l^? = A«[p«(e.)pW(efe)]-'/'exp[-(£W - £«)V(2An] (46) 
where Aq"^ and are constants and where 

p«(e)«p(''exp(A£) (47) 

is the level density and l3i — 1/Ti the inverse temperature of the heat bath 
labelled (z). 

The sum Vl^(') is also a Gaussian random variable with mean value zero. 
This fact suffices to derive for the average density matrix p{t,t') of the total 
system governed by H the integral equation 



p{t,t') = U{t)p{Q,0)W{t)+ I (It dr' U{t~T)Wp{T,T')W W{t' -t') . 

Jo Jo 

(48) 

The averaged time-evolution operator obeys 

U{t) = exp[-iHot]~ J dti J dt2 exp[-iHo{t - ti)] 

xW expi^iHoih - t2)]W Jjlpi) . (49) 
With the help of the assumptions (|45|) it is easily seen that 



{mla2b \ U{t) \ nlc2d) 

= 5™„5.Adexp[-z(i?™ +41) +ep))t- (l/2)(rW +r(^^)i]e(t) (50) 
where Q{t) is the Heaviside function and where 

r^l = 27r^(mza|W^W|m6)2 6{E„, + - E,, ~ 4'^) , z = 1, 2 . (51) 

nb 

In Eq. (|50p the sum of two Fs appears because W^^'^ and VF*-^-* are uncorrelated 
and may cause different damping. 

We take the trace of Eq. P5|) with respect to both heat baths and obtain 

^(mla2fe I p{t,t') I nla26) = ^cxp[-i(£;„< - E„t')] 

ah ah 

xexpH(4'^+4'^)(<-t')] 

X exp[-(i/2)(r« + T^^,)t - (i/2)(rw + r|f,))t'] 

x(TOla26 I p(0,0) I nla26) 

+ <^mnV I dT [ dT'expH(£;,„+£W+4'))(i-T)] 

Jo Jo 



ab 
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xexph(l/2)(rW +ri^))(t-r)] 

X (^\{mla2b I WW \ nlc2b)\'^{nlc2b \ p{t, t') \ nlc2b) 

^ nc 

+ ^|(TOla2& I VF(2) I nla2d)\^{nla2d \ p{t, t') \ nla2dU 

nd ' 

xexp[z(i?™ + £«+4'))(i'-T')] 

xexp[-(l/2)(rW +r|^))(t'-r')] . (52) 

The Kronecker delta on the right-hand side of Eq. (|52p follows from the first 
of Eqs. (pS)) . Eq. ([5^ shows that the nondiagonal elements of the reduced 
density matrix of system S decay exponentially in time. Therefore, we focus 
attention on the diagonal elements with 111 = 12, for t = t' denoted by Pn{t)- 
We differentiate Eq. ((52|) with respect to t and t'. In the resulting gain terms 
we use the weak-coupling assumption and the resulting Markov approximation. 
We make use of Eqs. ^ and (gT]). With 

W^liT,) = 2n4^\Q(:l\' exp[-(i?„ - i?„)V(2A?)] exp[(l/2)ft (£;,„ - E^)] 

(53) 

that yields Eqs. (P) and 
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